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Abstract 

In this study the magnetization phenomenon has been investigated 
as a behavior of interacting elementary moments ensemble, with the 
help of Ising model [1] in the frame of non-extensive statistical me- 
chanics. To investigate the physical systems with three states and two 
order parameters, the spin-1 single lattice Ising model or three states 
systems are used. In the manner of this model thermodynamical prop- 
erties of a great deal of physical phenomena such as ferromagnetism 
in bilateral alloys, liquid mixtures, liquid-crystal mixtures, freezing, 
magnetic orderliness, phase transformations, semi-stable and unstable 
states, ordered and disordered transitions [2,3,4,5]. 
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1 Introduction 



In this study, taking spin-1 Ising systems as the model dilute magnetic sys- 
tems have been investigated in the frame of non-extensive statistics [6,7,8,9]. 
The general Hamiltonian of next nearest pair interaction spin-1 Ising sys- 
tems is in the form 

H=-jJ2S i S j -Kj2sfS]- 

<ij> <ij> 
N N 

j=l j=l <ij> 

M^SiS^Si + Sj) (1) 

<ij> 

where J, K, D, H, L, M and N are bilinear exchange interaction constant, bi- 
quadratic exchange interaction constant, crystal field interaction constant, 
magnetic field due to S', dipole-quadropole interaction constant, magnetic 
perturbation of third degree and number of lattice points respectively. It 
the systems stay in semi-stable state or phase their properties change con- 
siderably. For examples, some alloys or metals lead to amorphons (glassy 
metal) structure when they are rapidly cooled down. This structure namely 
the phase is a semi-stable state. In this manner, some of the characteris- 
tics of the alloys or metals such as magnetic properties, resistance against 
corrosion, exhausting, wearing and hardness improve considerably. This 
Hamiltonian includes all possible terms when sf = Si it does not include 
higher order powers of spin. 

2 Dilute Systems 

Dilute magnetic system is a system which is obtained by inclusion of non- 
magnetic atoms to magnetic atoms. Exhibition of magnetic property of such 
systems is possible when the concentration of the magnetic atoms reaches 
to a certain value. Dilute magnetism is given by the expression A V B\- V C 
where A shows magnetic, B non-magnetic and C purely magnetic states. 
The Hamiltonian of such a system is defined in the form 

H= — ^2 SiSjJniUj (2) 

<ij> 
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which is obtained by taking M=K=L=0 in Eq.(l). In this equation Si 
and Sj are the spin vectors of i th and j th atoms, ni and rij are spin disorder 
variables taking the values and 1 whose average value gives the magnetic 
concentration c. J is exchange energy. 



Table 1: Spin Shaping Probabilities. 
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In table 1 where the spin shaping probabilities are presented, Vi are 
different shping numbers having same probability, en, ei2, £13, £21, £22, 
£23) £31; £32 and £33 are atomic interactions. If magnetic atoms are grouped 
as A and B denotes the atoms which are not magnetic then en = £13 = 
£33 = £AA,£i2 = £32 = £21 = £23 = €ABande 2 2 = cbb- The spin directions 
of the magnetic atoms are 1 and 3 and the directions of the atoms which 
are not magnetic is 2. The probabilities of the atoms having spins in the 
upward directions is x±, ,n downward is X3 and the atoms that are not 
magnetic is %2- Since the probabilities having symmetry are equal as a 
result y\2 = 2/21, 2/23 = 2/32 and 2/13 = 7/31 . The shaping of pairs at the lattice 
point takes place in nine different form and these shaping have been given 
in Tabic 1. 

3 Free Energy of the System 

The energy per atom of such a system is 
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= Jl{yi3 - VU + 2 V33 ) + ^AA{yu + 2yi3 + 2/33) + 2eAB + 

(Z/12 + 2/23) + tBBV22\ + 

(2/11 + 2/12 - 2/33 - ym)PH (3) 

where 7 is the number of nearest lattice points, N is the number of lattice 
points of the system, H is the external magnetic field and [3 = jr^f • Using 
the weight factor xi as the internal variable which has been developed for 
three state systems and y\j as double variable, the statistical weight of the 
spin-1 Ising system could be written as; 

w *= ^MT (4) 

Ll* Hn^=i(y^)!] 2 

where N is the number of lattice points in the system and 1 is the number 
of systems in the ensemble. On the other hand the definition of entropy is 

When the entropy is calculated for a single system (L=l) 
S 3 3 

— = fc B [( 7 - 1) Xilnxi - (|) Vij ln Vij] ( 6 ) 

i=l i,j=l 

is found. Thus free energy becomes: 



F = E-TS. 

Using the expressions in Table 1, free energy takes the form 



Tp T / 2/H + 2/33 X . 

F = J7(2/i3 g ) + 

7 

2 [ e ^A(yn + 2yi 3 + 2/33) + 2eAs(2/i2 + 2/23) + £#52/22] - 

1 3 7 3 

(2/11 + 2/12 - 2/33 - Vast)PH - - [(7 - 1) Xilnxi - ( -) Vijlnyij] (7) 

^ i=i «,j=i 
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4 Obtaining the Equilibrium State of the Dilute 

System 



The equilibrium state of a dilute system could be obtained by the calculus 
of variations. For this purpose three independent variables rj, £1 and £2 are 
defined: 

yi3 + 2/31 , s 

V = o = f 13 ( 8 ) 



2& = x 1 -x 2 (9) 



26 = J/12 - 2/32 (10) 

where 77 indicates the energy variation, £1 represents magnetization and 
£2 shows the abundance of the upward spins with respect to the downward 
spins of A atoms which are neighbors of B atoms. 



5 Equilibrium State of the Dilute System in the 
Frame of Non-extensive Statistical Mechanics 

Until 1998 all of the physical quantities of statistical systems were obtained 
by Boltzmann-Gibbs statistics. According to Boltzmann-Gibbs statistics, 
macro quantities such as free energy, entropy and internal energy of a sta- 
tistical system were accepted as extensive quantities. In 1998 a generaliza- 
tion in the Thermodynamical meaning has been carried aut to understand 
the structure or to solve a great number of unfamiliar systems and the 
generalization was inspired from the probability definition of the multifrac- 
tal geometry. Magnetization is a process with a long range and memory. 
Such systems are non-extensive. Thus magnetization will be investigated in 
the frame of non-extensive statistical mechanics. This generalization is the 
parametrization of all of the statistical quantities by a parameter q. In the 
limit q — ► 1 the statistics under investigation reduces to Boltzmann-Gibbs 
statistics. For the values of q different than 1, macro quantities such as 
internal energy, free energy and entropy in Boltzmann-Gibbs statistics are 
not extensive quantities in other words they are non-extensive. For gener- 
alization, mathematical expressions obtained from non-extensive statistical 
mechanics are introduced; 
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x l ~ q - 1 

ln„ x = and 

1-q 

exp* = {[1 + (1 - q)x]^ifl + (1 - q)x > 

Oelse. (11) 

On the other hand, when Eq.(lla) is written down for any function f(x), 
its derivative becomes 

1-q 1 ~ 1-q ~ f(x)l [U) 

In this system, the equilibrium state is determined using the variation 
of the free energy with respect to parameters 77, £1 and £2 given above. The 
variation of free energy leads to nonlinear equations: 

2.1 

2/112/33 = yi3 e q T ( 13 ) 



(2/ii)i = ( ^)7-i e p e (7-i)(^-^) 

2/33 ^3 

M = (^)2 e 2 (2/i 2 " < '-S/32" 9 )+(2/33" 9 -2/n" 9 ) ( 15 ) 
2/33 2/32 9 

Starting from Table 1 and taking x^~ q = x\~ 9 in Eq.(13) and 2/12" 9 — vl^ 9 
and y\^ q = y\^ q in Eq.(14) these equations are 
brought to simpler form: 

2J 

2/112/33 = 2/l3 e <7 ( 16 ) 

2/3 ff 

= (^i)7-i e ^ (17) 

2/33 £3 

Vll = [2/12^ (lg) 
2/33 2/32 

It is obvious that the solutions of these equations could be obtained 
depending on the entropy parameter q. In order to proceed to physically 
measurable results, let us relate the probabilities with the relevant concen- 
trations. Existence probability of the magnetic atoms in the lattice point is 
riA and the probability for the atoms that are not magnetic is ub- Therefore 
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n A = x 1 + x 3 orn A = y n + y 12 + y 13 + 2/31 + 2/32 + 2/33 (19) 

n B = x 2 orn B = y 21 + y 22 + y 3 2 (20) 

n AA = Vll + 2/13 + 2/31 + 2/33 (21) 

n AB = 2/12 + 2/13 + 2/32 (22) 

«BA = 2/21 + 2/23 (23) 

"-SB = 2/22- (24) 

Writing c = ija in other wards by taking the concentration of the mag- 
netic atoms above equations could be expressed in terms of c: 

n AA = c 2 (25) 

n AB = c(l - c) (26) 

n BA = c(l - c) (27) 

n BB = (l-c) 2 . (28) 

6 Calculation of the Orientation Probabilities 

In order to solve the equations given above a transformation given by 

- = ej* (29) 

is used in Eqs.(15),(16) and (17) which leads to the solutions of these 
nonlinear equations: 



6 = ^tanh 9 ( 7 -l)t. (30) 



The third parameter 



7 



£ 2 = ^tanh g ( 7 -l)i (31) 
is obtained. Using these results the probabilities X{ and yij are found: 

2(7-1)* 

2/n = —rr ( 32 ) 

2[coshg2( 7 - l)t + eg sT ] 

yi2=2/2i = ^(tanh ? ( 7 -l)t+l) (33) 

yi3 = 2/31 = —r ( 34 ) 

2[ex P g sT cosh, 2( 7 - l)t + 1] 

2/22 = n BB (35) 

2/23 = 2/32 = -^(1 - tanhg( 7 - l)t) (36) 

-j I 

2/33 = )t + exp 9 s . (37) 

2 expq sT [coshg 2( 7 — 1 
On the other hand the internal variables determined as, 

n A (l +tanh 9 7 t) 
37 = - (38) 

x 2 = n A (39) 
23 = reA(1 - t 2 anh " 7t) (40) 

7 Determination of the Physical Quantities 

When the relation between the parameter t and temperature is determined 
by using the obtained probabilities, the nonlinear equations and expression 
related to the concentration one gets: 

p vr n A tanh g7 t 



coshg 2(7 — l)t[nAA tanhq 2(7 — l)t + uab tanh g (7 — l)t — tia tanh 5 7^] 

(41) 



n J 4fitanh (? (7 - l)t 



coshg 2(7 — l)t[riAA tanhg 2(7 — l)i + uab tanh ? (7 — l)t — tia tanh g jt] 
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Magnetization 
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Figure 1: Variation of magnetization with respect to temperature for simple 
cubic structure (7 = 6.) at a certain concentration (c=l) and for different 
g values (Eq(29)). r = f(T) plot. 



According to Eq.(29) the magnetization is 

T q = 2^\orT q = n A tanh(qjt) (42) 

Variation of the magnetization with respect to temperature at a certain 
concentration value (c=l) for different q values, in other words the plot of 
the function r q = f(T) is given in Fig.(l). In the figure solid line refers to 
q = l, closer spaced dashed line shows the plot for q=0.8 and wider spaced 
dashed line for q = 0.5 . 

The susceptibility is expressed as 



j 

J_ _ kT[l + e$ T coshg 2(7 - l)t](n A - in AA - n AA ) + 2n AA (j - 1) 
Xq /? 2 [1 + ef cosh g 2(7 - l)t\{n\ + n AA n A ) - n A n AA 

In Fig. (2) the variation of susceptibility with respect to temperature at 
a certain concentration value (c=l) for different q values, in other words 
the plot of the function \q = f(T) is represented. In this plot solid line 
correspond to q = 1 and the dashed line to q = 0.5. 

When the magnetic part of the energy is taken under consideration one 
writes: 

^L = ll(4r ] -n AA ) (44) 



E q _ —~fJ 1 — ekT sec(foq)2(7 — l)t 
77 = ~ n AA[- — 



q _ _^ nAA[ - ~ : ^y» qJ -y, (45) 

1 + ekT sec(/i g )2(7 - l)t 



Using these expressions the specific heat is found to be; 
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Susceptibility 




Figure 2: Variation of susceptibility heat with respect to temperature for 
simple cubic structure (7 = 6).) at a certain concentration (c = 1) and for 
different g values (Eq(42)). x = f( T ) P lot - 
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Figure 3: Variation of specific heat with respect to temperature for simple 
cubic structure (7 = 6.) at a certain concentration (c=l) and for different 
g values (Eq(44)). SI = f(T) plot. 



_Cg_ = ^dE^ , } 

k B N k B N dT> v ' 

£± = l£ d ^ sec{h q )2(j - l)t - I 

kN 4k nAA dT [ e ^ sec{hq)2{l _ 1)t + l l [ > 

In Fig. (3) the variation of the specific heat with respect to temperature 
at a certain value of concentration (c=l) for different g values, in other 
words c q = f(T) plot is shown; where the solid line is for q = 1, closer 
spaced dashed line for q = 0.8 and wider spaced dashed line for q = 0.5. 



8 Conclusions 

In this study; by considering the dilute magnetic systems as an ensemble 
of interacting elementary moments and in this context with the help of 
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spin-1 Ising model; they are investigated microscopic level. Due to the 
statistical mechanics a bridge has been constructed between microscopic 
level approach and macroscopic experimental results. Magnetization is a 
long range phenomenon with a memory. Such systems are non-extensive. 
For this reason; magnetization has been considered in the frame of non- 
extensive statistical mechanics. Starting with the standard approach the 
generalization process proceeds. Same of the experimental studies where 
this type of systems are investigated has also been given [10,11,12]. In 
this study the variations of magnetization, susceptibility and specific heat 
with respect to temperature have been investigated. It is observed that; at 
different values of the entropic index q, in other words when q is decreasing, 
magnetization exhibits a linear variation rather than a parabolic one. On 
the other hand, it seen that the susceptibility does not undergo any change 
with q values. In the specific heat however, with decreasing q temperature 
dependence increases. 
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